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¢ Concept Schematic of tube-on-tank concept 
¢ Objectives of the CFD analysis 
¢ Details of the CFD model 
¢ CFD model results: 
¢ Incoming gas with temperature: 
¢ 273 K (warm case, baseline) 
¢ 100 K (cold case) 
¢ Tank wall boundary condition: 
a) Define wall temperature at 90 K 
b) Define the heat flux as 12 W/m2 


¢ Conclusions 
¢ Forward works 
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Introduction 


¢ The in-situ production of propellants for Mars missions 
will utilize Mars atmospheric carbon dioxide (CO.) to 
produce oxygen. 
¢ The oxygen is then cooled, liquefied, and stored to be 
available for Mars ascent propulsion system, which could 
be up to 2 years after liquefaction starts. 
¢ Recent investigations have demonstrated the feasibility 
of using high-efficiency reverse turbo-Brayton-cycle 
cryocoolers to: 1 stage tank 
¢ Cool the oxygen gas s | 
- Liquefy the oxygen gas sie ited 
¢ Achieve zero boil-off wy u 
¢ Control the pressure of S/4 5 . 
oxygen within a tank Helium tank < 2nd stage 
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A configuration of tube-on-tank liquefaction using a cryocooler. 
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¢ The gaseous Neon circulating in the cryocooler system is 
maintained slightly below liquid oxygen saturation temperature 
and is routed through a network of cooling tubes. 

¢ The oxygen gas produced from the in-situ production process is 


introduced into the chilled tank. 
7/25/2016 TFAWS 2016 — August 1-5, 2016 5 


eds . Objectives of the CFD analysis NASA 


¢ Model the proposed active thermal liquefier design 
(tube-on-tank) to reduce the uncertainty of the heat 
transfer coefficient inside the tank 
¢ Model liquetier configuration options to create an 
efficient system 
¢ Meet the requirement of the liquefaction time 
¢ Minimize the mass and power for the active 
thermal liquefier system 
¢ Understand the relationship between the incoming 
gaseous O, temperature versus tank surface area 
and condensation rate of the gas inside the tank 
¢ Investigate the advantage of pre-chilling gaseous O, 
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Jr FD model approach using ANSYS Fluent NASA 


¢ 2D axisymmetric, nc = 11201, nv 
11785, dt=0.01s 
o Solver: Pressure-based, 
transient, coupled 
o Multiphase model: 
Mixture/slip Velocity/Implicit 
Body force 
o Turbulence model: shear 
stress transport (SST) k-w (2 see 
eqns) Eee 
¢ Tank wall boundary condition ¢ : 
a) Set temperature: 90 K 
b) Set heat flux: 243.6 
W/20.3 m2 = 12 W/m? 
(based on the lift of the 
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¢ Inlet gaseous O.: 
¢ Warm gas at 273 K (baseline) 
¢ Cool gas at 100 K 
¢ Mass flow rate: 
¢ 2.2 kg/hr = 6.11 g/s (baseline) 
¢ Initial conditions: 
¢ T=100K 
¢ VOF of the liquid oxygen: 0 
¢ Inlet tube: 
¢ 1” diameter at the top of the tank 
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4 ANSYS Fluent results: t = 38 mins, case (a) NASA 


Temperature 
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Contour 
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¢ Liquefaction occurs at the bottom of the tank 
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¢ Free convection inside the tank and near the 
interface of liquid and gaseous O, 
¢ Flow streamline contours shown 


7/25/2016 TFAWS 2016 — August 1-5, 2016 10 


/XANSYS Fluent results: t = 38 mins, case NASA 
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¢ Heat transfer coefficient near the dry wall is around 50 W/m2-K. 

¢ Natural convection calculated from CFD model is an order of 
magnitude (50 W/m?-K) larger than hand calculations using 
Grashof numbers (0.3 — 1.85 W/m?-K) 
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Bede ANSYS Fluent results 


Time history of the volume-average gas temperature 
o Incoming gas: 273 K and 100 K 
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¢ In graph (a), with tank wall at 90 K, GOX chills down very quickly, 
within 10 minutes for both cases - incoming gas at 273 K and 100 K. 
This is the optimal case. 
¢ In graph (b), the tank wall heat flux is fixed. This is the worse case. 
¢ With the incoming gas of 273 K, it takes much longer to cool the 
gas down and the gas is much warmer. 


¢ With the incoming gas of 100 K, it takes 20 minutes to chill down. 
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Time history of the mass of Lox 
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¢ Graph (a) wall temperature fixed at 90K; graph (b) heat flux is fixed. 
¢ The liquid oxygen inside the tank at t = 40 minutes Is 
¢ For incoming gas of 273 K: 
¢ 1.48 kg in case (a), 0.55 kg in case (b), a factor of 2.7. 
¢ For incoming gas of 100 K: 
¢ 1.52 kg in case (a), 0.95 in case (b), a factor of 1.6. 
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ANSYS Fluent results 


Temperature contour of mixture 
o Incoming gas: 273 K 
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¢ Graph (a) wall temperature fixed at 90 K; graph (b) heat flux is fixed. 
¢ The warm gaseous O, chills down within smaller volume with a cold 
wall. 
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ANSYS Fluent results 


Temperature contour of the mixture 
o Incoming gas: 273 K 


9.20e+0 siesta 

9.19e+0 exer? 

9.18e+0 OMSex0) 

9.17e+0 (a) O:1fe*0 

9.16e+0 9.168+0 (b) 

9.158+0 9.15e+0 

9.14e+01 9.146+0 

9.13e+0 9.13e+0 

9.12e+0 9.126+0 

9.10e+0 9.10e+01 

9.09e+0 9.09e+0 

9.08e+0 9.08e+01 

9.07e+01 9.07e+0 

9.06e+01 9.06e+0 

9.05e+0 9.05e+01 

9.04e+01 9.04e+0 

9.03e+01 9.03e+0 

9.02e+01 9.02e+0 

9.016+01 9.016+0 

9.00e+01 9.006+01 ¥ 
8.99e+01 8.99e+0 

Contours of Static Temperature (mixture) (k) nie Senaee 03) Jun 14, 2016 Contours of Static Temperature (mixture) (k) (Time=2.4000e+03) 5, 2016 
Fluent Release 16.2 (axi, dp, pbns, mixture, sstkw, transient) ANSYS Fluent Release 16.2 (axi, dp, pbns, mixture, othe “iat sient) 


¢ Graph (a) wall temperature fixed at 90 K; graph (b) heat flux is fixed. 
¢ With a cold wall, condensation occurs more surface area of the tank. 
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ANSYS Fluent results 


Volume fraction of liquid contour 
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¢ Graph (a) wall temperature fixed at 90 K; graph (b) heat flux is fixed. 
¢ With a cold wall, more liquid O, is formed at the bottom of the tank. 
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\ Conclusions NASA 


¢ The condensation rate predicted by CFD analysis depends on the tank 
wall boundary condition. We presented two bounding cases. 


— For case (a), with an assumed 90K tank wall the gas chills down very quickly, within 10 
minutes for both incoming gas of 273 K and 100 K. This is the optimal case. 


— For case (b), with a constant heat flux assumed, for the incoming gas of 273 K, the 
condensation rate is much smaller and larger portion of the tank area near the top stays 
warm. This is could be the worse case. 

— Heat transfer around the majority of tank is natural convection driven. 

— Incoming warm gas induces mixing currents and forced convection occurs near the inlet 
tube. 

¢ Natural convection calculated from CFD model is an order of magnitude 
(50 W/m?-K) larger than hand calculations (0.3 — 1.85 W/m?-k) 


¢ Tube-on-tank concept works for the baseline condition (warm case). MAV 
tank provides enough heat transfer area for liquefaction. There are some 
concerns near the top of the tank, however, it is beyond 95% fill level. 

¢ Pre-chilling gaseous O, will speed up the liquefaction rate inside the tank 
as long as the lift of the cryocooler is allocated enough for the tank itself. 
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¢ 1D thermal model in Matlab and 3D thermal model in Thermal 
Desktop are on-going to include the cooling fluid in the model 
to get a more realistic tank wall boundary conditions. 


¢ Test plans on the tube-on-tank concept are on-going. The 
objectives of the test plan are 
— To integrate the reverse turbo-Brayton-cycle cryocoolers system with 
the tank and control system (existing hardware) 
— To Investigate the performance of the active thermal liquefier system 
under 
¢ Different fill levels 
¢ Different feeding gaseous oxygen temperatures 
¢ Different control schemes 
¢ Validations between the test results and model results will be 


performed. 
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